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In chaotic reaction-diffusion systems with two degrees of freedom, the modes governing the 
' exponential relaxation to the thermodynamic equilibrium present a fractal structure which can be 

characterized by a Hausdorff dimension. For long wavelength modes, this dimension is related to the 
Lyapunov exponent and to a reactive diffusion coefficient. This relationship is tested numerically on a 
reactive multibaker model and on a two-dimensional periodic reactive Lorentz gas. The agreement 
with the theory is excellent. KEY WORDS:Reaction-diffusion systems; hydrodynamic modes; 
fractals; microscopic chaos. 
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I. INTRODUCTION 
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C/3 ^ Recently, many studies have been devoted to fractal structures induced by the chaotic dynamics in the phase-space 
of non-equihbrium statistical systems and, in particular, to the fractal character of the hydrodynamic modes governing 
the exponential relaxation of the system to the thermodynamic equilibrium [0-|2j . This fractal character has been 
related to the entropy production in the approach of equilibrium, for a diffusive multibaker model Q . Previous works 
' had shown similar results for the entropy production in non-equilibrium steady states On the other hand, the 
(— ( , fractal dimension characterizing the hydrodynamic modes of diffusion has been related to the diffusion coefficient 
Q ■ and to the positive Lyapunov exponent by a relation, first obtained for diffusive multibaker maps and proved for 
O ' general chaotic systems with two degrees of freedom [Q. The construction of the hydrodynamic modes in Refs. 

gives us a new approach to relate macroscopic transport coefficients and microscopic chaotic properties. It establishes 
' a link between the irreversible transport processes and the reversible microscopic dynamics. 
^ , Two other approaches are the so-called escape rate formalism and the thermostated system approach. The first 
' approach can be applied to open systems with absorbing boundaries [^-^. The escape of trajectories leads to the 
I formation of a fractal repeller, consisting in the set of orbits forever trapped within the absorbing boundaries. Under 
such conditions, the transport properties can be related to the positive Lyapunov exponents and the Kolmogorov-Sinai 
' entropy per unit time as well as to fractal dimensions |^-^. In the second approach, an external field is applied to the 
system, and in order to keep the kinetic energy constant, a special force acts on the particles as a heat pump [|o|-^. 
The system is still time-reversible but does not preserve phase-space volumes anymore. The trajectories converge 
"Y"; 1 to a fractal attractor. In thermostated systems, the transport properties have been related to the sum of Lyapunov 

■ exponents ||ll|,|l^. In the new approach of Refs. the advantage is that the fractal curve considered is directly 
^ , related to the hydrodynamic modes of relaxation towards thermodynamic equilibrium. It does not require absorbing 

I ■ boundaries neither a thermostat. 
■ O ' In the present paper, we extend the new approach of Refs. to reaction-diffusion processes. Reaction-diffusion 
^ ■ processes are of fundamental importance to understand self-organization in physico-chemical systems ||l^ , [l^ . In 
^ ' chemical systems, the microscopic mechanisms of the relaxation toward the thermodynamic equilibrium are still 
poorly understood already for simple reactions. For the purpose of contributing to this important question, we 
consider here a simple reactive process of isomerization. The particle is supposed to perform a deterministic motion 
among static scatterers and to carry a color A or B. When colliding on some special scatterers, playing the role of 

■ catalysts, its color changes instantaneously with a given probability pa, 

a 

A + catalytic scatterer < — > B + catalytic scatterer , with probability po , 
A + other scatterer < — > A + other scatterer , 

B + other scatterer < — > B + other scatterer . (1) 

For this class of systems, we obtain an expression of the Hausdorff dimension of the reactive modes of relaxation in 
terms of the reactive dispersion relation and of a function Q{a,f3) generahzing the Ruelle topological pressure P(/3). 
In the long wavelength limit, we then infer a relation between the Hausdorff dimension of the modes, the reactive 
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diffusion coefficient, the reaction rate, and two derivatives of the function Q{a,P). In the hmit pq — > 0, we recover 
a relation previously derived for the diffusive case [^[jj]. This new expression relates the macroscopic transport and 
reaction processes to the microscopic underlying dynamics. 

The paper is organized as follows. In Section^, we generalize the work of Refs. |^,^ to the reactive case considered 
here. In Section III, we test our relation on a reactive multibaker model. We next study the case of a two-dimensional 
periodic reactive Lorentz gas in Section IV . Conclusions are drawn in Section M. 



II. RELAXATION MODES: GENERAL RESULTS 



We consider a system such as a Lorentz gas in which the particles move independently of each other along trajectories 
of a deterministic dynamical system given by 

The variables of this system are for instance the positions and velocities of the particle: X = (r,v). The vector field 
F(X) is defined in this space. The system is spatially periodic in the positions r. The particle moves in a periodic 
lattice of scatterers given by a certain potential of interaction. 

Moreover, each particle carries a color which is either A or B. The color changes when the particle interacts with 
some of the scatterers, called the catalytic scatterers according to the reaction (|l|). We notice that we may also 
consider that the particle has a spin one-half which flips upon collision on special scatterers. 

We can model the reactive events by a change of color with probability po when the trajectory meets a certain 
hypersurface S(X) = defined in the space of the variables X. This hypersurface surrounds each catalytic scatterer 
and constitutes the locus where the reaction occurs. The introduction of the color has for consequence that the phase 
space is composed of two copies of the space of the variables X, which are glued together along the hypersurface E. 
This hypersurface can be viewed as the gate from copy corresponding to the color A to the one corresponding to color 
B and vice versa. These two copies now form two subspaces of the whole phase space of the reactive system. In this 
framework, a reactive event corresponds to the crossing of the hypersurface S from one color subspace to the other. 

A statistical ensemble of particles of both colors is described by two probability densities 5(X, A, t) and g(X, B, t), 
to find a particle at the point X with color c = ^ or _B, at the current time t. The time evolution of these probability 
densities is ruled by two coupled Liouville equations: 

dt (z(X, A, t) + div [F(X) g(X, A, t)] = -Po S [I](X)] [g(X, A, t) - g(X, B, t)] , , . 

dt g(X, B, t) + div [F(X) g(X, B, t)] = +po S [S(X)] [g(X, A, t) - g(X, B, t)] . 



where s is a special distribution defined as 



S[E(X)] /i(X) = lim 



F(X).g(X) 



«[S(X)] A|X-fF(X)l (4) 



6 being the Dirac distribution. The distribution (Q) expresses the fact that the particles change their color upon 
crossing of the hypersurface S. Accordingly, the probability densities just after the hypersurface E in the direction of 
the vector field F are the following linear combinations of the densities just before E: 

g(X+,A,0 = (l-po)'z(X_,A,i)+po9(X_,B,i) , 

g(x+,s,<)=pog(x_,At) + (i-Po)g(x_,B,<) , 

where X± = X* ±eF(X*) with E(X:„) = and an arbitrarily small e > 0. These linear combinations can be expressed 
in terms of the transmission matrix 

J=(^-P^ /° V (6) 
\ Po 1-poJ 

We denote by X = $*Xo the trajectory of the differential equation (||) from the initial condition Xq. Moreover, 
we gather both probability densities in the vector 
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The solution of the coupled Liouville equations (^) from the initial probability densities qo can then be written in the 
form of the following Frobenius-Perron operator: 

/- \ T^t(X) . On f^^'Xl 



9X 



with the matrix (g) and where Nt (X) is the number of crossings through the hypersurface E performed by the segment 
of trajectory from Xq ~ $~*X to X = $*Xo. If the flow of Eq. (|^) is volume preserving, we have that divF = 



which implies that 



1, which is the case for the systems considered in the following. 
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The system of Liouville's equations (|3|) defines a random process and it is important to characterize the dynamical 
randomness generated by this process. For this purpose, we calculate its Kolmogorov-Sinai entropy per unit time which 
is the amount of information produced per unit time by the process. The Kolmogorov-Sinai entropy is calculated 
by partitioning the whole phase space into small cells and by computing the probability for a trajectory to visit 
successively different cells of the partition If we suppose that the differential equation (||) defines a chaotic 

dynamical system with Lyapunov exponents! A^}, the dynamical instability of the deterministic motion is a first 
contribution to the Kolmogorov-Sinai entropy given by the sum of positive Lyapunov exponents according to Pesin's 
theorem | [T5[ . A second contribution comes from the random change of color upon crossing of the hypersurface E. At 
each crossing, the particle jumps from one color subspace to the other with probability pq. If the number of crossings 
is equal to Nt during the lapse of time t, the frequency of S crossings is equal to : 

u ^ lim (9) 

t-+oo t 

where (•) denotes the average with respect to the equilibrium state. At each crossing, we have an independent dichoto- 
mous random variable with probabilities {pq, 1 — po), which defines a Bernoulli process with the rate v. Accordingly, 
the overall Kolmogorov-Sinai entropy per unit time of the process (|^) is given by 



Xi>0 



Poln— + (1 -po)ln- 



Po 1 - Po 



(10) 



Pesin's expression of the Kolmogorov-Sinai entropy is recovered in the deterministic limits po — and po = 1. 

For < Po < Ij we remark that the Kolmogorov-Sinai entropy ( p^ remains positive and finite as for a fully 
deterministic dynamical system. The way the reaction is modelled here is therefore compatible with a deterministic 
dynamics that would take place near each catalytic scatterer on a time scale much shorter than the free flight between 
each scatterer. We can imagine that the reaction follows a complicated trajectory on such a short time scale in 
a complicated potential, which is here simplified by capturing the complexity of the reaction into a sole transition 
probability Pq. Still the so-defined random process has a finite Kolmogorov-Sinai entropy contrary to stochastic 
processes which have an infinite Kolmogorov-Sinai entropy per unit time. Indeed, it has been shown elsewhere 
1^,0 that stochastic processes such as the birth- and-death processes or the linearized Boltzmann processes have 
a dynamical randomness on arbitrarily small scales e with an e-entropy per unit time increasing as h(e) ^ ln(l/e). 
Stochastic processes defined by Langevin stochastic differential equations have an even larger dynamical randomness 
with an e-entropy per unit time increasing as h{e) ~ The Kolmogorov-Sinai entropy is the supremum of the 

e-entropy in the limit e — > and is therefore infinite. In contrast with such stochastic processes, the present random 
process has a finite Kolmogorov-Sinai entropy. This dynamical randomness is comparable with the randomness of 



the classical dynamics on graphs which has been recently characterized |18 . For the classical graphs, the motion is 
one-dimensional so that the only source of randomness is the random bifurcations at the vertices of the graph. Here, 
there is an extra source of randomness which is the dynamical instability due to the chaotic motion of the particle in 
the multi-dimensional phase space of Eg. (^) . 

The overall dynamical randomness induces a relaxation of the probability densities toward the thermodynamic 
equilibrium. In order to study this relaxation we have to find the relaxation rates which are given by the eigenvalues 
of the Frobenius-Perron operator (||). It turns out that the Liouville equations can be decoupled if we introduce 
the probability density / to find a particle at the point X no matter which is its color, and the difference g between 
both color densities: 

f{X,t)^q{X,A,t)+q{X,B,t) (11) 
g{X,t)^q{X,A,t)-q{X,B,t). (12) 

These new functions are gathered in the vector 
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qd(X,t)= ( {(J'fl ) (13) 



,3(X,t) 

which is related to the previous vector ^ by the relation 

Qd = S • q (14) 

in terms of the matrix 



^1 -1 

This matrix diagonalizes the transition matrix m) into 



Td = S-T.S-i=(i , ). (16) 



l-2po 



Accordingly, the Liouville equations (||) decouple into 



ft/(X,i)+div [F(X) f{X,t)]=Q, 

dt g{X, t) + div [F(X) .g(X, t)] = -2po S [E(X)] ^(X, t) 



(17) 



which shows that the probability density / obeys the deterministic Liouville equation while the difference g solely 
controls the reaction. Similarly, the Frobenius-Perron operator (H) splits in the two following operators: 

/(x,t)^ (18) 



g(x,t)^(i-2po)^^(^ \--:r^ T' . (19) 



^JVt(X)_ffO_(*_j2^ 

l^(*-*x) 

An eigenvalue problem can now be considered separately for both operators (|8|) and ([T9|). If we denote them 
generically by P*, we are looking for eigenstates given by some distributions which are quasiperiodic in the position 
space r with a wavenumber k: 

*k(X) = e'''"-(^) Tk(X) (20) 

where Tk is periodic in r. The wavenumber k describes spatial inhomogeneities of wavenumber i = 27r/||k||. The 
distribution (pQ) is supposed to be a solution of the eigenvalue problem 



P*^'k = e"''*^'k (21) 

associated with the operator P*. The quantity — Sk is the relaxation rate of the eigenstate (po|). 

The idea is that these eigenstates can be obtained in principle by applying the operator P* over a long time to 
some expression (GQ) and by renormalizing the result in order to compensate the relaxation according to 



4'k(X) = lim e^""-' P' 

t^OG 



^^k.r(X) Y(x) 



(22) 



even starting from some arbitrary function T. This renormalization-semigroup idea has been developed in detail for 
diffusion in Refs. ||l^,U|. For chaotic dynamical systems (||), such eigenstates turn out to be singular distributions 
without a density function but with cumulative functions defined by 

Pk(^) = / d(?'*k [X(0')] (23) 
Jo 

where 6* is a parameter (such as an angle) of a curve X(6') defined in the phase space. 

The hydrodynamic modes of diffusion as well as the reactive modes can be identified as the eigenstates corresponding 
to the slowest relaxation rates and they can be constructed thanks to the above idea of renormalization semigroup, 
as we shall show below. 
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A. Diffusive modes 



In order to construct a diffusive mode of wavenumber k we consider tlie Frobenius-Perron operator ( [l8| ) and an 
initial distribution T(X) of particles of initial positions Tq. We integrate the trajectories by using the equations 
of motion (H) in order to obtain the position after a long time t. As a consequence of the aforementioned idea 
of renormalization semigroup, the relaxation rate of the diffusive mode is given by the decay rate of the Van Hove 
intermediate incoherent scattering function as 

Sk = lim 7 ln(exp[ik • (rt - ro)]> (24) 



where (•) is an average over initial conditions distributed according to some distribution T p] |l7|p4[ |. In a similar way, 
the cumulative function of the diffusive mode can be obtained as 

, „„, J^O-exp(.l.-|,.(0-)-r„(.-)|) 
/j" df exp{tk ^ |r,(e') - r„(»0|} 

where the initial conditions are distributed according to a distribution T which is here taken as a uniform distribution 
along a one-dimensional line parametrized by the angle G [0, 27r[ 

For two-degree-of-freedom chaotic systems with one positive Lyapunov exponent A, the cumulative function ( [25| ) 
forms a fractal curve in the complex plane and its Hausdorff dimension is given in terms of the diffusion coefficient 
and the Lyapunov exponent according to 

L>H(k) = 1 + ^fc2 + 0(^4) (26) 
A 

as shown in Refs. The derivation of Rcf. Q is here summarized because we need it for the generalization to the 

reactive modes. 

We assume that we can use the result by Sinai, Bowen and Ruelle that averages can be performed in terms of 
unstable trajectories covering the phase space within a certain resolution pl|-p3|. The probability weight given to 

each of these trajectories is inversely proportional to their stretching factor Aj which characterizes their dynamical 
instability. Accordingly, Eq. (|24|) can be transformed into the condition 



^ \Ai'^\-' exp(-ski) exp{ik • [rj^'^ - v^'^]} ^t^^ 1 . (27) 



j 



If the sum in Eq. ( [2^ ) is restricted to the trajectories issued from initial conditions in the interval [0,9], we obtain 

at time t a polygonal approximation of the cumulative function (p5|): the integral is an average over trajectories 
with initial conditions in [0,9] and the denominator is proportional to the factor exp(ski). At time t, the curve 
(Re .Fk, Im fk) C C is thus approximated by a polygon of sides given by the small complex vectors 

AF(^') = |ApV exp(-Ski) exp {ik • [rp'^ - r[,^'^] } . (28) 

Each side has the length 

= |AF(J)| = |A|^V^ exp(-Re Sk t) . (29) 



In the limit t ~f oo, this polygon converges to a fractal curve, characterized by a Hausdorff dimension given by 
e^^ ^ 1. Accordingly, the Hausdorff dimension of the hydrodynamic mode should satisfy the condition 

J2 lApV'' eM-Dn Re t) ^t^^ 1 . (30) 
j 

On the other hand, Ruelle's topological pressure is defined in dynamical systems theory by 

P(/3) = lim i H[At\'-^) (31) 

t^oo t 

where the average (•) is carried out over the equilibrium invariant measure p3t . The mean positive Lyapunov exponent 
of the system is given by A = —dP/df3[fj=i and in this system without escape, P{1) = 0. Equation (|3^) can be 
transformed into 
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^lAp'r'' exp[-P(/3)i]^,^o,l. (32) 
j 

Comparing with Eq. (p7|), we obtain 

P{Dh) = Dh Re sk . (33) 

If the wavenumber k vanishes, the cumulative function ( |25| ) becomes i^k=o(^') = ^'/27r, which forms a straight Hue 
in the complex plane. In this equilibrium limit, the relaxation rate vanishes, Sk=o = Oj so that P{Dii) = 0. It 
means for this system without escape that Dn = 1, as it should be for a straight line. For a non- vanishing but 
small wavenumber, the Hausdorff dimension is expected to deviate from unity. Inserting Z^h — 1 + S and the 
dispersion relation Sk = — + 0{k^) in Eq. (|33|), we can expand both sides in powers of the wavenumber, using 
the aforementioned properties of the topological pressure. This straightforward calculation shows that the Hausdorff 
dimension of the hydrodynamic mode is given by Eq. (^6|) . 



B. Reactive modes 



1. General results 

We are now going to extend this result to the reactive modes of the isomerization process (|l|). The modes of 
interest here are the ones describing the evolution of the difference between the two colors A and B, a quantity which 
relax to zero when the system reaches equilibrium. In this case, we consider the operator ( p^ ) which differs from the 
diffusive Frobenius-Perron operator ( [T^ ) by the extra factor (1 — 2po)^*- Accordingly, this factor should multiply the 
exponential exp(ik • r) in the expression ( p^ in order to obtain the relaxation rate of the reactive modes, which is 
thus given by 

5k = lim I ln((l - 2^0)"^' exp[ik • (r^ - Vq)]) (34) 

where Nt is the number of collisions performed on catalysts during the time t ||26|| . The cumulative function of the 
reactive modes are similarly given by |^ 

, 1,„ />»'(l-2.o)".<'-'e,M.k. [...(.•) -r„(r)l} 
-~ dO' (1 - 2p„)".l«') exp{jk ■ [r,(90 - r„(e')]} 



We notice that Eqs. ( |34| ) and (|35|) consistently reduce to Eqs. (^4|) and ( |25| ) in the purely diffusive limit po 0. 
As for the diffusive case, Eq. (|34|) can be transformed into the condition 

^ |A«|-i (1 - 2po)^"' exp(-Skt) exp{*k • [rj^") - r'^")]} 1 . (36) 

j 

The curve (Re _fk, lui -Fk) C C is approximated at time i by a polygon of sides given by 

e, = |AF(^)| = lAp'^l-i |1 - 2por"' exp(-Re Sk t) . (37) 



In the limit t — > oo, this polygon converges to a fractal curve, characterized by a Hausdorff dimension given by 
J2j ^ 1- Accordingly, the Hausdorff dimension of the reactive mode should satisfy the condition 

J2 lAp^r^" |1 - 2por"'^« exp(-i?H Re Sk t) ^t^^ 1 . (38) 

3 

By analogy with Eq. (|3^), we can define a function Q(a,/?) by 

J2 lAp^r'^ |1 - 2^01""^"' exp [-Q{a,m ^^^^ 1 (39) 

3 

which can be rewritten as 
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Q{a,f3) = lim 1 ln((l - 2^0)""^' \At\'~^} ■ (40) 



Comparing Eq. (39) with Eq. (|3^), we obtain 

Q{Dii,DH) = DiiResk. (41) 
The relation between this new function and the topological pressure is given by 

P(/3) = g(0,/3). (42) 
If the wavenumber k vanishes, the cumulative function (|3^) is reduced to 

F^=o{0) = hm J.. ...... ■ (43) 



This is not a straight line as in the diffusive case: each initial condition is here weighted by a factor (1 — 2po)^*- 
Provided that po < 0.5, this means that the contribution of a given initial condition to Fq will be smaller if the 
corresponding trajectory meets many catalysts, and vice versa. But the Hausdorff dimension of _Fo is anyway equal 
to unity. Indeed, for k = 0, the relaxation rate is given by 

Sk=o = lim i ln((l-2po)^*), (44) 

t^oo t 

which is also equal to (3(1,1) according to the definition (^0|). Using these equalities, we can deduce from ( |4l| ) 
that Dh = 1 for vanishing wavenumber. Again the Hausdorff dimension will deviate from unity for a non-vanishing 
wavenumber. The reactive dispersion relation is given by = — we compare it with ( ^ ) at 

vanishing wavenumber, we obtain 

Q(l, 1) = -2k = lim I ln((l - 2^0)"^') ■ (45) 

Using this equality, the dispersion relation and the expression D}i{k) = 1 + ak^ + 0(fc''), Eq. ( plj ) can be expanded 
in powers of the wavenumber, giving an expression of the Hausdorff dimension of the reactive mode 

^H(k) = 1 - ,^^og.f'^.^%9.i + ■ (46) 



An expression for the two partial derivatives of Q{a,(3) can be directly obtained from Eq. (|40|), provided that 
Po < 0.5: 



and 



^^n^)- lim ((l-2po)^'l^|A«l) f,^. 



Similarly, the reactive diffusion coefficient is given by 



((1 - 2po)^* (Vt - To)') 

2dt({1^2po)^*) 

where d is the dimension of the space of positions r where diffusion occurs 



t^oo 2 d t {{1 ^ 2pa)^*) ^ ' 
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2. Limit of small reaction probability po <^ 1 



Let us now consider the limit case po ^ 1- The three terms of the denominator appearing 
expanded in powers of po- The reaction rate is given by 



2k = hm - 

t — *oo t 



2po {Nt)+2pl {Nt) + -pl {Nt) 



-2pl {{Nl)^{N,f)-^pl {{Nl)-{N,f) 



2K = 2vpo + 2{v-a^)pl + 4(- v - a'' +t)pI + 0{pt) 



in terms of the frequency (^) of coUisions on the catalysts v and of the quantities 



a-=^hm - {{N^)-{Nt)^) 



lim - (\{Nf) + \{^N,f-{^N^){N,) 

t— ►CO t \ ^ 6 



The two further terms are given by 



and 



where 



lim — 

t — >oo t 



2po {Nt)+2pl {Nt) + -pl (Nt) 



^Apl {(Nf) (N,)') ^8pl {(Nf) (N,)^) 
+ 12 pI (^l{N!) + l{N,r~{N^) (N,) 

-2k + 2a^pl + 4{a^ ~ 2 t) pi + 0{p^) 



|f(l,l) = -A + 2po lim Ic^xit) 

op t^oo t 



+ 2 pI hm - [(1 + 2{Nt)) C^xit) - D^x{t)] 



t^oo t 



4po lim - 

t^oo t 



- + 2{Nt) + 2 {Nt)^ - {N^] 



-{l + {Nt)) D^xit) + -E^xit) 



Oipt) 



dQ 

9/3 



(1, 1) = -A + a^xpo + P^xpI + 7«aPo + C(Po) 



A= hm - (In I A, I) 

t^oo t 

C, x{t)^{Nt ln|At|)-(7Vt)(ln|At|) 

D, x{t)^{N! ln|A,|)-(iV,2)(ln|A,|) 

E, x{t) = {Nf ln|A,|)-(iV3)(ln|A,|) 

a„A = 2 lim - Ci^x{t) 

t^oo t 

P^x = 2 lim l[{l + 2{Nt))C^x{t)~D^x{t)] 



7«;A = 4 hm - 



+ 2(Nt) + 2 (Nt)^ - (Nf) C^x(t) 



0{pt) 



C^x{t) 



-{l + {Nt))D^x{t) + -E,x{t) 



For these coefficients, the index k refers to the reaction, related to the number of collisions on a catalyst Nt, while 
the index A refers to the dynamical instability, described by the stretching factor At and the Lyapunov exponent A. 
The reactive diffusion coefficient (MS) can also be expanded in powers of the reaction probability po as 



=I? + Apo + 0(Po) 



(62) 



in terms of the standard diffusion coefficient 



V = hm 



1 



t^oo 2 d t 



((r.-ro)^) 



(63) 



and of the coefficient 



A = — lim — 

t^oo d t 



[{N,{r,-ror)-{N,){ir,-rof)] 



(64) 



which characterizes the correlation between the transport by diffusion and the reaction. 

We notice that, in the limit po 0, the reaction rate (50) vanishes while the reactive diffusion coefficient 
to the diffusion coefficient (|63|), as expected. 

Using (p3[) and (p3), the Hausdorff dimension (E6|) becomes 



tends 



i^H(k) - 1 



~ Po 



- ^ Po 



2 r 



Po 



a „ 
2-pl 



A2 



Po 



A3 



A2 



"^pI 



(To To 



pI + 0{p%) 



+ o(fc*) . 



(65) 



or 



^H(k) 



1 



O(po) 



A:2 + 0(fc4) . 



(66) 



We immediately see that for po ^ 0, Eq. ( |65| ) reduces to the expression ( p6[ ) obtained for the diffusive case. 

We are now going to test this formula for two examples: a reactive triadic multibaker and a reactive two-dimensional 
periodic Lorentz gas. 



III. THE REACTIVE TRIADIC MULTIBAKER MAP 



A. Description of the model 

The multibaker map is a simple model used to mimic a diffusive motion P,|25|] . In the case of the triadic symmetric 
multibaker, we consider an infinite one-dimcnsional chain of squares, 

Si^{{x,y,l): 0<x<f, 0<2/<f}, (67) 

where Z e Z is the label of the square. A point particle performs jumps from square to square, according to the 
transition rule Z (/ — 1) if a; < 1/3, Z — > Z if 1/3 < x < 2/3 and Z (/ + 1) if a: > 2/3. In order to introduce a 
reaction, the point particle carries a color A or i3, and the central part |] of one square over L is catalytic: if the 
point particle is in this portion of the square, it changes its color with a probability po (see Fig. The map of the 
model is thus given by 

r (3a;, f,/- l,c) , 0<a;<l/3 
<p(x,y,;,c) = i (3a;- 1,2^,;,£(0 c) , 1/3 < a; < 2/3 (68) 
i (3a;-2, 2^,; + l,c) , 2/3 < a; < 1 

where c is a discrete variable taking the values +1 or —1 when the color of the point particle is respectively ^ or _B, 
and 



9 



^/^x _ f — 1 if Hs a catalytic square, with a probability Pq , /ggx 
I +1 otherwise . 

We consider the probability density q{x, y, Z, c) which evolves under the Frobenius-Perron equation 

(qt{^,3y,l + l,c) , 0<y<l/3 
qt+,{x,y,l,c)^ I [l-po{l)]qt{^,3y-l,l,c)+po{l)qtm^,^y-lJ,-c) , 1/3 < j/ < 2/3 (70) 
[qtm^,3y~2,l~l,c) , 2/3<t/< 1 

where 

/>\ _ ( po if / is a catalytic square , 
P°^''~\0 otherwise. ^ ' 

The Kolmogorov-Sinai entropy per iteration (|l^) of the reactive triadic multibaker is given by 

/^KS=ln3 + ff Inl + i-I^ln-^ (72) 
3L Po 3L 1 - Pq 

because the frequency (^ is equal to the ratio of the area of a reactive rectangle to the area of L squares: v = 1/3L. 
We notice that the result (|7|) can also be obtained by using a generating partition which establishes the equivalence 
between the reactive multibaker map ( |68|) and a Markov chain (see Ref. |]l7|). 
If we define the difference between the densities of the colors A and B, 

gt{x,y,l) = qt{x,y,l,+l) - qt{x,y,l,-l) , (73) 

we obtain the following equation for its evolution 

r 32/, / + 1) , 0<2;<l/3 

gt+i{x,y,l)^l[l-2po{l)]gt{'^,3y-l,l), 1/3 < y < 2/3 (74) 
Ut(^, 32/ -2,^-1) , 2/3<y< 1 

which is the form of the reactive operator ( |l9| ) in the case of the present multibaker model. 



B. Relaxation rates and eigenmodes of the reactive operator 



We restrict our description to a unit cell of length L, I = being the catalytic square, and we choose quasiperiodic 
boundary conditions, assuming that the solution of the equation ( [74[ ) is quasiperiodic on the chain with a wavenumber 
k. Moreover, we suppose exponentially decaying solutions, with decaying factor x(fc) = exps^ where |x(^)l < 1 or 
Resk < 



gt{x,y,l) - X* exp{ikl). 
The following operator Qh governs the time evolution of gt{x, y, I) p 



Qk 



' gt+i{x,y,0) 
gt+i{x,y, 1) 

9t+i(x,y,L- 2) 
gt+i{x,y,L- 1) 



= - 2/)5t(f , 32/, 1) + (1 - 2po)eiy - i)f?(| - y)gt{^,3y - 1, 0) 

+e-''=^% - |)5t(^, 32/ - 2, L - 1) 

= - y)5t(f , 32/, 2) + % - i)0(| - y)gt{^, 3y - 1, 1) 

+^(y-|)3t(^,3j/-2,0) 



= 0(1 - y)gt{l,3y, L-l) + 0{y~ i)0(| - y)gtm^,3y -l,L-2) 
-%-f)3t(^,3j/-2,L-3) 

. -e»'=if?( 1 - 2/)5t(f , 32/, 0) + 0iy - i)0(| - y)gtmi,3y - 1, L - 1) 



f).gt(^,3j/-2,L-2). 



We now want to obtain the eigenvalues and the eigenstates of this operator 



(75) 



(76) 



(77) 
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We suppose that the leading eigenstates are uniform along the unstable direction x 
If we now cumulate over x and y, we get 



where 



Gt{x,y,l)= dx' dy'^{x',y',l) = xC{y,l) 
Jo Jo 



C{y,l)= / dy'y{y',l) 



(78) 



(79) 



(80) 

is the reactive cumulative function (|35| ) for the multibaker. We replace gt{x, y, I) by y{y, I) and gt+i{x, y, I) by x3^(2/i 
in Eq. {u&j and we integrate over the interval [0,?;]. We obtain 



C(y,0) = 



Ciy,l) 



C(3y, 1) 



< 2/ < 1/3 

[C(l, 1) + (1 - 2po)C{3y ~ 1, 0)] , 1/3 < y < 2/3 

4 [C(l, 1) + (1 - 2po)C(l, 0) + e-'^^^Ciiy - 2, L - 1)] , 2/3 < y < 1 



^C'(3y,2), 

[C(l,2) + C(3y-l,l)] 



< 2/ < 1/3 
1/3 < y < 2/3 



\- [C(l, 2) + C(l, 1) + G(iy - 2, 0)] , 2/3 < y < 1 



3x 



C(y,L-2) 



3t^C(3y,L-l), 

[C(l,L-l) + C(3y-l,L-2)] 



< y < 1/3 
1/3 <y< 2/3 



\- [C(l,L-l) + C(l,L-2) + C(32/-2,L-3)] , 2/3 < y < 1 



3x 

^C^(3y,0), 



C(y, L - 1) = <( 3L [e^fcL ^(10) ^ c[iy - 1, i - 1)] , 



< y < 1/3 
1/3 < y < 2/3 



}- [e^'=^ C(l,0) + C(l,L-l) + C(3y-2,L-2)] , 2/3 < y < 1 



3x 



The eigenvalues are obtained by setting y = 1 in Eq. (81), which leads to the eigenvalue equation 



/(l-2po-3x) 1 
1 (1 - 3x) 1 

1(1- 3x) 



V 





Akh 







(1 - 3x) 
1 



ikL 



\ ( ^^(1,0) \ 
C(l,l) 
C(l,2) 



. C(l,i-2) 

(i-3x)y Vc(i,i-i)y 



= 



(81) 



(82) 



We thus have to solve the determinant of an Hermitian matrix to find x(^)- For H^is purpose, we assume a solution 
of the form 



C(l,/) = Ae+''^'+Be-*^' 
After substitution in the system (|8^), we find that 

1 + 2cos6' 
X = 7, 



(83) 



(84) 



and 



sm6'L 

cos kL = cos vL + po ^ ■ 

smfc' 



(85) 



For the present reactive system, the relaxation rate s(\t) = Inx(fc) depends on the wavenumber fc as s(fc) = —2k 
_^0(fc4), where K is the reaction rate and I?^"^^ is the reactive diffusion coefficient. 



11 



Solving Eq. ( p5| ) for its leading root, we can thus obtain the reaction rate and the reactive diffusion coefficient either 
in powers of the transition probability pq: 

I>"' = 5 + ^ + 0«), (87) 

or in the asymptotic limit L — > oo: 

In the limit po ^ 0, Eq. (^6|) shows that the reaction rate vanishes linearly with the transition probability and Eq. 
( ^ ) that the reactive diffusion coefficient tends to the purely diffusive one V = 1/3, as expected. On the other hand, 
Eq. (|8^) shows that the reaction rate decreases as when the catalytic sides are diluted in the limit L ^ oo, as 
expected for a diffusion-controled reaction in a one-dimensional system Moreover, the vanishing of the reactive 
diffusion coefhcient ( ^9| ) as for L ^ oo is consistent with the results previously obtained in Ref. We notice 
the remarkable property that, in the diffusion-controled regime of the reaction, the leading term of the asymptotic 
expression for the reaction rate ( ^8|) does not depend on the transition probability po, which only appears in the next 
terms. Therefore, the leading term of ( p8| ) no longer expresses the property that the reaction rate vanishes in the 
limit Po — > 0, as shown by Eq. (|8^). This is consistent with the fact that the reaction rate is largely independent of 
the local reaction rate at the catalysts in the diffusion-controled regime, as already noticed in Refs. |p|,p6[. 
Let us consider more explicitly the simple case L = 1. In this case, Eq. (^2|) reduces to 



which gives the solution 



For fc ^ 1, x(fc) is then given by 



3x - 2 cos fc - 1 + 2po = (90) 



f,. l-2po + 2cosfc 
X{k) = ^ ■ (91) 



X(fc) ~ e-^-e-^'"'' 

^Xo(l-2?^^^fc'). (92) 



Expanding ( |9l] ) in powers of k up to the second degree, we obtain 



l-TT^^ (93) 
3 - 2po / 

which allows us to identify xo = = e^'^" and P^') = 

The dependence of the reactive diffusion coefficient P'''^ on po is shown in Fig. ^ for different values of the chain 
length L. The value of the diffusion coefficient V is indicated for comparison. 



C. Fractal dimension of the cumulative functions of the eigenmodes 

We now want to study the cumulative functions C{y,l), solutions of Eq. (^l|), and more specifically to calculate 
their dimension. These functions are represented in the complex plane in the case L = 1 in Figs. ^ and ^. Fig. ^ 
shows their dependence on the reaction probability Pq, the wavenumber being fixed, k — 0.1. Their dependence on 
the wavenumber for po = 0-3 is shown in Fig. ^. As explained in Section ||, these complex cumulative functions can 
be approximated by a polygon at a given time t, which means at a certain iteration step n in this discrete case. This 
concept will here be applied using a symbolic dynamics. 

The symbolic dynamics is constructed by partitioning the phase space of the multibaker maps depicted in Fig. 
Each of the L squares of the chain is divided in three parts for x € [0, a; |] or x 1]. These 3L rectangles 
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are successively assigned the symbols ujn G {0, 1, ■ ■ ■ ,3L — 2,3L — 1} from left to right along the chain. At the n 



th 



step of the construction described in Subsection II B , the cumulative function is a polygon in the complex plane with 
its sides corresponding to a subrectangle associated with a sequence of n symbols {ti;„_iaj„_2 ■ • -^^i^a}- Moreover, 
we notice that the sequence of symbols also determines the square I where the point is located along the chain. At 
the next steps of the construction, the cumulative function keeps its values already obtained at the borders of each of 
the aforementioned subrectangles. We can thus denote the values of the cumulative function at the lower and upper 
borders y±(ujn^i, L0n^2, ' ' ' , wi, loq) of one of the subrectangles in the square I more shortly by 

C±(w„_i, w„_2, . . . Wo) = C [y±(w„_i,a;„_2, . . . , wi,a;o), I] ■ (94) 
The complex vector defining the side of the polygon at the n*'' step of the construction is thus given by 

AC(W„_1,CJ„_2, ■ • ■ Wo) = C+(cJ„_i, W„_2, ■ • . , Wi, Wo) - C_(w„_i, W„_2, • • • , Wi, Wo). (95) 

The succession of two symbols wiwo corresponds to a transition in the multibaker phase-space, with which a certain 



probability is associated . Using 



we can construct a 3L x 3L transition matrix M of elements Afi^^i^g describing 



the factor connecting the portion of square corresponding to wo to the one corresponding to wi. Using this matrix, 
we can describe the time evolution of the intervals AC(w„_i, w„_2, • ■ • , wi, wq) as 



AC(Wn_i, W„_2, . . . ,Wi,Wo) = 



■ A/^2i^j Af^j^o AC(wo) 



(96) 



We define 



E 



^1, W„_2, . . . , Wi, Wo)| 

E (-^""')c.„-,c.o |AC(wo)|^ 



D 



|AC(wo)|^ 



(97) 
(98) 



where the matrix M. is composed of the elements: 

Supposing that this matrix admit the spectral decomposition, A4 = /i (//|, Eq. (| 

r^(C) = EH^) (a^H 



(99) 



can be rewritten as 



(100) 



where u denotes the vector of unit elements — 1 and v the vector with the elements v^^ = |AC(w)p. According 



to the definition of the Hausdorff dimension, for n — > oo, the quantity (100) must converge to a finite non-vanishing 
value if _D = Dn. This will only be the case if = 1. Z?h is thus solution of the equation obtained by imposing 1 as 
eigenvalue of the matrix A4. Let us consider explicitly the case L = 1. We have to calculate the 3x3 determinant 



- 1 









1 


1 


3x 




3x 









- 1 




l-2po 


l-2po 


l-2po 


3x 




3x 


3x 













1 


1 


1 


3x 




3x 




3x 



= 0. 



(101) 



This gives the equation 



3X 



1 = 0. 



(102) 



We will solve (102) in the limit case po 

jjiO) _ In 3 



0, using the expression of x in the case L ^ 1 (|91|). For po = 0, we have 
H ~ in(i+2cosfc) • other hand, the Hausdorff dimension takes the unit value for k = 0. Substituting the 

expansion Dh = 1-1- ak^ + 0{k'^) in Eq. (102) and expanding in powers of the wavenumber fc, we can solve for the 
unknown coefficient a and obtain 
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1 



(3 - 2po) ln(3 - 2po) - (1 - 2po) ln(l - 2po) 
Now expanding in powers of the transition probability po, we finally obtain 



InS 3 



9^ 



27 



1 



(In 3)2 



4 2 
27^0 



(103) 



(104) 



We can again compare this expression with the general result ( |65|) of Section For the triadic multibaker, the 
Lyapunov exponent is given by A = In 3. Moreover, the stretching factor In |Af | is constant in the triadic multibaker, 
which implies that the coefficients C„a, I?kA: -E'kA and therefore a^A, PkXi luX defined in Eqs. (|5^)-(|6l]) are vanishing. 
Whereupon, Eq. (65) is here reduced to 



D^{k) = 1 



A 



1 



3 



Using the expansion of the reactive diffusion coefficient up to the third order in po 

1 4 



Eq. (|105|) becomes 



-DH(fc) = 1 



1 1 
hTs 3 
1 1 



1 



3^° 



9 



27 



(105) 



(106) 



(In 3)2 3 



- 4 fc^ - 



(In 3) 



27 
1 
3 



(In 3)2 
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2 ^8rp3 fc2 + 0(p^A:2) + 0(fc4) 



(107) 



Let us now 



The second term of the right members of (104) and (107) are equal. They correspond to ^ 
compare the other corrections in po- The only correction in in (107) is jy^^ ^ 2 Pq k^. We can evaluate 



= hnif 



{{Nf) — (iVt)2) in the case of the triadic reactive multibaker L — 1. Indeed here, colliding on a 



catalyst corresponds to being in the region x gIj, |] 
then 



The mean number of collisions on a catalyst per unit time is 



1 



(Nt) 







and 



which finally gives us 



1 

3 

1 

3 



1' + 



(108) 



(109) 



(110) 



The coefficient of the correction in pg is then equal to (^i^-^^y^ ^ , 
correction in p^ in (107) contains three different contributions: 
Again we can evaluate r 

1 /I 



which is exactly the coefficient obtained in (104). The 



(In 3)^ 9 



p3 fc2 + 



(InS)^ 3 



p3 fc2 



(hTsp 



r pI k"^ 



= hm - - {Nt) + -{N,Y~{Nt){Nt) 



1 

27 



2 

81 



(111) 



We obtain 

Eq, 



80 



^^^g^ -j J which is the coefficient expected for the correction in p^ . Accordingly, we have confirmed that 

^ |) for the Hausdorff dimension of the reactive modes is valid for the reactive multibaker map. The validity of 

the approach is furthermore tested in Fig. ^, where the Hausdorff diemansion _Dh is plotted versus the square of th e 
wavenumber k'^ . The symbols are obtained numerically by solving Eq. (102) and compared with the expression (104). 
The agreement is good up to po = 0.3. We notice that, beyond this value, further terms would be required in the 
Taylor expansion (104). 



14 



IV. THE 2D REACTIVE PERIODIC LORENTZ GAS 



We are now going to consider as second example a two-dimensional reactive periodic Lorentz gas. This model has 
been introduced in Ref. Q. In this system, a point particle, carrying a color A or B undergoes elastic collisions 
on hard disks fixed in the plane and forming a regular triangular lattice characterized by the lattice fundamental 
vectors ei = d{l, 0) and 62 = d{^, ^). The radius of the disks is assumed to be equal to unity and d is the distance 
between their centers. We shall work in the finite horizon regime, 2 < d < for which the diffusion coefficient is 
finite. Some of the disks play the role of catalysts: when the point particle collides on one of them, it changes its 
color instantaneously with a probability po- The catalysts form a regular triangular superlattice over the disks lattice, 

with as fundamental vectors Ei — nd{0, — \/3) and E2 — nd{^, ^), where n is an integer parameter controlling the 
density of catalysts: in the direction of Ei and E2, one disk over n is a catalyst. Globally, one disk over N — is 
a catalyst. The configurations with n — 1 and n = 2 are depicted in Fig. ^. 

Our goal here is to study the cumulative reactive eigenfunctions of this system. As done in Ref. j26j, we will 
construct the Frobenius-Perron operator governing the time evolution of the color density describing the reactive 
process in this system. We will then perform a spectral analysis of this operator. 



A. The reactive operator 



The variables needed to describe the point particle in this system are its position, its velocity and its color, i.e. 
{x,y,Vx,Vy,c), where c is a discrete variable taking the values +1 or —1 when the particle color is A or B. However, 
the collisions on the disks being elastic, the energy of the point particle is conserved. The magnitude of its velocity 
is a constant of motion, that we suppose equal to one, ||v|| = 1. The velocity can therefore be described by only one 
coordinate, which can be for example the angle between the velocity and the x-axis. This system is periodic, which 
means that instead of describing the dynamics in the full space, we can restrict our description to an elementary cell 
of the superlattice of catalysts, containing N disks (see Fig. Moreover, the relevant information concerning the 
dynamics is contained in the collision dynamics: between two collisions, the point particle performs a free flight. The 
three-dimensional flow dynamics can therefore be reduced to a Poincare-Birkhoff map, describing the dynamics from 
collision to collision. We will use the Birkhoff coordinates x = (j, 6*, tu), 1 < j < iV = 3n^, < ^' < 27r, — 1 < < 1, 
where j is the index of the disk of the elementary cell on which the collision takes place, 6 is an angle giving the 
position of the impact on this disk and zu is the sine of the angle between the outgoing velocity and the normal at 
the impact. In these coordinates, the mapping of the reactive Lorentz gas is given by 

x„+i = v(x„) , 

— ~t~" ^(^n) ; ('112') 
In+l = In + a(x„) , 
Cn+1 = £(X„) C„ . 

tn is the time of the n"^ collision. T(x) is the first-return time function. The vector 1„ gives the cell of the superlattice 
visited at time i„ and a(x„) is a vector giving the jump carried out on the superlattice during the free flight from x„ 
to x„+i. The change of color is controlled by the function: 

/ N _ r —1 if jn is a catalyst, with a probability po , (US) 
" I -f 1 otherwise . 

We notice that the dynamical randomness of the reactive Lorentz gas is characterized by the following Kolmogorov- 
Sinai entropy per unit time 



h-KS = X + I' 



Po In — -I- (1 - Po) In ■ 



Pq 1 - Po 



(114) 



where A is the unique positive Lyapunov exponent of the two-dimensional Lorentz gas and ly is the frequency of 
encounters with a catalytic disk. 

We now want to obtain the Frobenius-Perron operator associated with this map. We first start from a description 
of the particle in the full space. The current state and position of the particle are given by a suspended flow in terms 
of the coordinates (x,r, 1, c), where < r < r(x) is the time elapsed since the last collision, 1 is a vector of the 
superlattice and c — ±1. The phase-space probability density depending on these variables g(x, t, 1, c) evolves in time 
under the Frobenius-Perron operator defined by Eq. (m pfl]. This continuous- time operator can be reduced to an 
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operator describing the dynamics from collision to collision under the mapping (112). The details are presented in 
Ref. . The operator for the modes of wavenumber k is given by 



Rks 9k,s(x, c) = exphsr(vp-ix) - ik ■ a(vp^ix)] (115) 
X {[1 -p(y)"^x)](7k,s(¥'"^x,c) + p(¥5"^x) (7k,s(¥'"^x,-c)} 

where 

p{^^{j,9,m})^\P" ifj is a catalyst, ^^^g^ 
L otherwise. 

The operator Rk^s is similar to a Frobenius-Perron operator of the mapping <p up to important factors which take 
into account the varying time of flight between the collisions, the spatial modulation by the wavenumber k, as well as 
the color change on a part of the domain of definition of the Poincare-Birkhoff mapping. If we now consider the total 
phase-space density / (x) = q{x, A) + q{x, B) and the difference of phase-space densities g{x.) — q{x, A) ~ q{x, B) , their 
time evolutions are decoupled and we can write an operator governing the evolution of (7(x), describing the reactive 
modes of wavenumber k in this system [ p6[ 

VFk,. ff(x) = [1 - 2p(y-ix)] exp[-sr((p-ix) - ^k • a{^-^x)] gi^-^^) , (117) 
which we call the reactive operator. 



B. Spectral analysis of the reactive operator 



To identify the eigenmode of the reactive operator (117) which controls the slowest relaxation of the color, we 
consider the eigenvalue problem for W^k,Sk and its adjoint imposing 1 as eigenvalue 

VFicsk V'k(x) = Vk(x) , (118) 

wis, V^k(x) = Vik(x) , (119) 

assuming the normalization condition 

(^kV'k) - 1 . (120) 

A formal solution for T/^k (respectively -(/jk) can be obtained by applying successively Wk,s (respectively M^k^) to the 
unit function: 

n 

M^) = lim n [1 - 2p(<p-™x)] exphskT(cp-™x) - ik ■ a(v-™x)] , (121) 

m— 1 
n-1 

V;k(x)= lim ]T [1 - 2p(<^^"x)] cxp[-4T(<p^^^x) + zk • a((^^"x)] . (122) 
The normalization condition gives an equation to obtain S]^: 

n 

1 = hm ( IT [1 - 2p(¥>-"x)] exp[-SkT(¥'-™x) - ik ■ a(v:""x)]) . (123) 

7n——n-\-l 

We notice that if we consider segments of trajectories of time t from initial conditions Vq we have that 

n 

^ T(c^-™x), (124) 

m— — n+1 
n 

Tt-Toc^- a(¥'"'"x) , (125) 

m— — n+1 
n 

(1 - 2po)^* ^ n [1 - 2p(<^"'"x)] , (126) 



m— — n+1 
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so that Eq. (123) can be transformed into Eq. (^J) for the reactfon rate. 

We are interested in the smallest relaxation rate — Sk which dominates at long times t. For the reactive case, the 
dispersion relation is given by 



Sk = -2k - I?('')k2 + ©(k^^) , (127) 

where k is the reaction rate and I?*^'-* is the reactive diffusion coefficient. Solving Eq. (|l23| ) allows us to obtain the 
values of these two coefficients. Fig. ^ shows the dependence of the reactive diffusion coefficient V^^^ on the reaction 
probability po, for two different configurations of catalysts, where n = 1 and n = 2. The value of the diffusion 
coefficient T> is indicated for comparison. We notice the remarkable result that the behavior is similar to the one 
observed for the multibaker model (see Fig ^). The explanation is the following. The case n = 1, with a single 
catalytic disk in each cell where the motion is ballistic, is similar to the multibaker map for i = 1 (see Fig. 0). In 
this case, the reactive diffusive coefficient rapidly increases with the transition probability pq in both the reactive 
multibaker model with L — 1 and the reactive Lorcntz gas with n = 1. In contrast, the reaction starts to be controled 
by the diffusion as soon as > 2 in the reactive Lorentz gas [ p6[ and, similarly, for the reactive multibaker with 
L > 2. In these cases, the reactive diffusion coefficient has a milder dependence on the transition probability pq and 
becomes even smaller than the diffusion coefficient for large values of po- 



C. Fractal dimension of the cumulative eigenmodes 



The cumulative function of the reactive eigenmode is obtained by integration of ( |121[) over the angle 
initial position and the a:-axis 



Fi,{0;j, n7) = — / de' lim [] [1 - 2p(v-™x')] exp[-Sk r((p-"x') - ^k • ai^- 

m=l 



between the 



(128) 



which is equivalent to the formula ( p5[ ) according to the correspondences (124)-(I26). 

Fig. 1^ shows these functions in the complex plane and as a function of the wavenumber k, in the case n = 1, 
d = 2.3 and po = 0.3. The Hausdorff dimension of these curves can be computed by a box-counting algorithm for 
different values of the wavenumber k. The values obtained are plotted versus in Fig. 0, in the case n = \, d = 2.3, 
Po = 0.1 and pa = 0.3. We have numerically computed the coefficients a^x, /3kA and cr, using the expressions (|5S|), 
(|^) and (|l|). We have obtained a^x ^ 0.23, ~ 0.18 and P^x ^ 0.15. The sofid fines in Fig. | are given by Eq. 
(pS), considering the corrections up to the second order in pq. We observe the remarkable agreement between the 
theoretical prediction and the numerical results. 



V. CONCLUSIONS 



In this paper, we have studied in detail simple reaction-diffusion systems in which an isomerization is induced by 
chaotic dynamics. We have constructed the reactive modes of relaxation toward the thermodynamic equilibrium. 
These modes are given by singular distributions having a cumulative function at small values of the wavenumber. The 
cumulative functions form fractal curves in the complex plane because of the chaotic dynamics and in spite of the 
randomness introduced by the transition probability po of the color in the isomerization. We explain this remarkable 
result by the fact that the randomness of the system is of a mild type since its Kolmogorov-Sinai entropy per unit 
time ( p^ remains finite. 

We have here derived a formula given by Eq. ( p6| ) for the Hausdorff dimension of the reactive relaxation modes in 
chaotic reaction-diffusion systems with two degrees of freedom (|4^). We have defined a function Q(a,/3) generalizing 
the Ruelle topological pressure P{P). For small values of the wavenumber, we obtain a relationship between the 
Hausdorff dimension Djj, the reactive diffusion coefficient 'D^^\ the reaction rate k, and two derivatives of Q{a,P), 
namely ^(1, 1) and ^(1, 1). In the limit of a small reaction probability po, the Hausdorff dimension is related to 

I?^'-* and to the Lyapunov exponent A, with corrections in powers of po- For po — 0, we recover the relationship 
( p6| ) obtained in Ref. for the diffusive hydrodynamic modes. We have shown that the reaction rate, the reactive 
diffusion coefficient and the Hausdorff dimension can be expressed in terms of quantities which characterize the 
statistical correlations between the number Nt of encounters with the catalysts, the distance Vt — Tq travelled in the 
diffusive motion across the system, and the stretching rates ln|Aj| characterizing the dynamical instability and the 
chaos. Our formula ( ^6| ) thus establishes a new relationship between macroscopic reaction-diffusion processes and the 
microscopic dynamical chaos. 
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We have tested our results on two reaction-diffusion systems with a simple isomerization reaction: a triadic multi- 
baker map and a two-dimensional reactive periodic Lorentz gas. In both cases, we study the limit of small wavenumbers 
and small reaction probabilities. The agreement with the theory is excellent. 

Moreover, the study confirms and clarifies results previously obtained for similar models of reaction-diffusion 
P, p5|j2^ . In particular, our results confirm that the reaction-diffusion process will obey on a large spatial scale 
and a long time scale the macroscopic reaction-diffusion equations: 

dp A „ d^pA , „ d^pB , . 

-^-T^AA-g^^VAB^^~<PA~PB). (129) 
dpB „ d'^PA , „ 9Vb , / ^ ,-,qn^ 

^ PgA-gp- -^n{pA-pB)- (130) 



for the color densities pA and pB- According to results (|62|) and (|64|), the diffusion coefficients of Eqs. (13C) are given 
by 

Vaa=Vbb = ^ =I?+-po + 0(Po) (131) 

V — P'"'^ A 

Vab=Vba = ^ = --po + 0(Po) (132) 

which explains that the coefhcients of cross-diffusion can be small if the transition probability po is small. Eq. ( |6^ ) 
shows that this cross-diffusion depends on the correlation between the number iVj of catalysts encountered and the 
square (rt — ro)^ of the distance travelled by the particle in its diffusive motion. Moreover, Figs. || and and Eqs. 
( p7| ) and (^) for the multibaker model show that the coefficient of cross-diffusion can become significant for high 
transition probability po and if the density of catalysts is very high, in the case L = \ for the multibaker and n = 1 
for the Lorentz gas. We may here say that the cross-diffusion, i.e. the coupling between the diffusions of colors A and 
-B, is induced by the reaction. 

In conclusion, the modes of relaxation toward equilibrium turn out to display unexpected fractal properties which 
have their origin in the chaotic dynamics of the microscopic motion of particles, even if this motion is time-reversal 
symmetric and volume-preserving. This result has been previously obtained for diffusive processes and, here, we 
have shown that this result extends to reaction-diffusion systems where it allows us to understand the relaxation to 
equilibrium at the level of the microscopic motion. 
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FIG. 1. Geometry of the reactive triadic multibaker map for various values of the distance L between the catalysts. 
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FIG. 4. Three-dimensional representation of the cumulative functions of the reactive eigenmodes of the triadic multibaker; 
(ReFfe,ImFfc) are represented in the complex plane as a function of the wavenumber k, in the case L = 1, po = 0.3. 
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FIG. 5. Reactive triadic multibaker: dependence of the HausdorfT dimension Dn of the cumulative functions of the reactive 
eigenmodes on the square of the wavenumbe r fc^, for L = 1 and for different values of the reaction probability po- The symbols 
are obtained by numerical resolution of Eq. (102) and the solid line is the analytic solution (104) obtained in the limit po <^ 1. 
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FIG. 8. 2D reactive periodic Lorentz gas. Three-dimensional representation of the cumulative functions of the reactive 
eigenmodes; (ReFk,ImFi^) are represented in the complex plane as a function of the wavenumber k, in the case n = 1, d = 2.3 
and po = 0.3. 
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FIG. 9. 2D reactive periodic Lorentz gas. Dependence of the Hausdorff dimension Dh of the cumulative functions of the 
reactive eigenmodes on the square of the wavenumber k'^, for n — 1, d = 2.3, po = 0.1 and po = 0.3. The symbols are obtained 
by a box-counting algorithm and the solid lines are given by Eq. (^5|) taking into account the corrections up to the second order 
in po. 
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